TO APPEAR IN " ClRCUMSTELLAR MEDIA AND LATE STAGES OF MASSIVE STELLAR EVOLUTION (2006)" RevMexAA(SC) 



SUPERNOVA EXPLOSIONS IN WINDS AND BUBBLES, WITH 
APPLICATIONS TO SN 1987A 
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RESUMEN 

Favor de proporcionar un resumen en espanol. If you cannot provide a Spanish abstract, the 
editors will do this. Massive stars can significantly modify the surrounding medium during their lifetime. 
When the stars explode as supernovae, the resulting shock wave expands within this modified medium and 
not within the interstellar medium. We explore the evolution of the medium around massive stars, and the 
expansion of the shock wave within this medium. We then apply these results to understanding the expansion 
of the shock wave in the ambient medium surrounding SN 1987A, and the evolution of the radio and X-ray 
emission in this case. 

ABSTRACT 

Massive stars can significantly modify the surrounding medium during their lifetime. When the stars explode 
as supernovae, the resulting shock wave expands within this modified medium and not within the interstellar 
medium. We explore the evolution of the medium around massive stars, and the expansion of the shock wave 
within this medium. We then apply these results to understanding the expansion of the shock wave in the 
ambient medium surrounding SN 1987A, and the evolution of the radio and X-ray emission in this case. 
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1. INTRODUCTION 



- SHOCK WAVES — ClRCUMSTELLAR MATTER 
SUPERNOVAE: INDIVIDUAL (SN 1987A) 



SUPER- 



Core-collapse supernovae (SNe) arise from mas- 
sive stars (> 8Mq). These stars tend to modify the 
medium around them substantially, often leading to 
the formation of circumstellar (CS) wind-blown bub- 
bles around the star. At the end of its life when 
the star explodes as a SN, the resulting shock wave 
will interact with the modified environment before it 
reaches the interstellar medium (ISM) which in some 
cases may take a substantially long time. The sub- 
sequent evolution of the shock wave, and the radia- 
tion signatures from the supernova remnant (SNR) 
will depend crucially on the structure and dynamics 
of the surrounding medium. In this paper we ex- 
amine the structure of the CS environment around 
stars, and the aftermath of stellar explosions which 
occur within this environment. Finally we apply the 
lessons learned to one of the most fascinating and 
well-studied objects, SN 1987A. In jg] we study the 
formation of the ambient medium around massive 
stars. In $3] we discuss the evolution of SNe in this 
medium. Finally in |2] we focus on applying these 
results to understanding the evolution and radiation 
signatures from SN 1987A. 
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2. THE ENVIRONMENTS OF MASSIVE STARS 

The environment of a massive star depends on 
its zero-age main-sequence mass, its rotational veloc- 
ity, the metallicity, and the presence of other nearby 
stars, among other factors. Therefore the evolution 
can be quite complicated, but several common fac- 
tors exist. There has been significant discussion of 
the evolution of massive stars by Norbert Langer at 
this conference, and a highly comprehensive review 
of observations of the circumstellar medium (CSM) 
around massive stars was given by You-Hua Chu. 
Therefore I will concentrate here on those details 
that are essential to understanding the evolution of 
the subsequent SN shock wave upon the explosion of 
the star. 

2.1. Main-Sequence Stars 

Massive stars generally start their lives as main- 
sequence O/B stars. A star of solar metallicity dur- 
ing its lifetime will usually lose mass through a ra- 
diatively driven wind (Kudritzki & Puis 2000) with 
a mass-loss rate on order 10 -8 to 10 _6 M© yr _1 , and 
with a wind velocity on order a few thousand km s _1 . 
The interaction of the wind from the star with the 
surrounding medium will lead to the formation of an 
interstellar wind-blown bubble (WBB), the structure 
of which was first delineated by Weaver et al. (1977, 
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hereafter W77). If the wind parameters are fixed, 
the bubble shows essentially 4 different regions as 
we go outwards in radius (1) Freely expanding wind. 
For a given mass-loss rate M and wind velocity v w 
the density of the wind goes as p w = M /(4irr 2 v w ). 
(2) A low-density, high-pressure (and therefore high 
temperature) region of shocked stellar wind. (3) A 
region of shocked ambient medium. The shocked am- 
bient medium usually cools quickly, resulting in the 
formation of a thin, dense shell. (4) The unshocked 
ambient medium. Most of the volume is usually oc- 
cupied by the region of shocked stellar wind, forming 
a hot, low-density cavity. W77 found a self-similar 
solution for the bubble with constant wind parame- 
ters, with radius varying with time as R cx t 3 / 5 . 

Several factors may complicate this simplistic de- 
scription. The wind properties are not constant but 
change with time as the star evolves through dif- 
ferent phases. Hydrodynamic instabilities and the 
onset of turbulence may cause a change in the dy- 
namical properties. Mixing at the interface between 
the hot shocked wind and the dense, cool, shocked 
ambient medium, which may also be a conductive 
interface, tends to lower the temperature in the in- 
terior. Density inhomogeneities or an asymmetric 
wind may lead to bubbles that are not spherical. 

The radius of the outer shock of the bubble (the 
outer edge of the dense shell) is (Weaver ct al. 1977): 

R. h = 0.76 (jj i 3 / 5 (1) 

where L = 0.5-MwJ i s the mechanical luminosity of 
the wind, M is the wind mass-loss rate and v w is the 
wind velocity. In the case of a main-sequence (MS) 
star with M = lCT 7 A/ yr" 1 , and velocity 2500 km 
s _1 , L is about 1.984 x 10 35 ergs s _1 . If this lasts 
for about 10 million years, the total energy released 
will be about 

E MS = 6.25 x 10 49 M_ 7 w 2 2 500 t 10 ergs (2) 

where M_7 is the mass-loss rate in terms of 1O -7 M0 
yr _1 , U2500 is the wind velocity in units of 2500 km 
s _1 and time is in units of 1 million years years (tio 
= 10 xlO 6 years). 

We assume that the main sequence star is formed 
in a medium with an average density of about 2.34 
x 10~ 23 g cm~ 3 (a number density ~ 10 parti- 
cles cm~ 3 , appropriate for an ionized region). From 
equation [1] the radius of the swept-up shell will be 

Rms = 48.8 M 1 ^ 5 v%l Pm 1/5 t 10 3/5 P c (3) 

where p is in units of 2.34 x 10~ 24 g cm~ 3 . How- 
ever we note that often, especially for less massive 



stars whose lifetime is large, the bubble will come 
into pressure equilibrium with its surroundings at an 
early stage, and thereafter will stall. The radius of 
the bubble may then be up to an order of magnitude 
smaller than is given in equation [3] (see Dwarkadas 
2006c, 2007). 

The ratio S TO of the mass swept up by the shell 
to the total mass in the wind is given by: 

_ / M ah \ _ 4n R 3 sh p a _ 4rr 0.76 3 f v w H 2 p \ 2 
m \M wmd ) 3 Mt 3 2 3 / 5 \ M ) 

(4) 

where p a is the ambient density. For the MS this 
becomes 

S m = 1.54 x 10 5 v% 5 00 t% 5 p 2 ^ MZ 2/5 (5) 

The wind velocity decreases slowly so the average 
velocity over the entire main sequence phase could 
be lower than that used. The duration of the main 
sequence phase also depends on the mass of the star, 
decreasing with increasing mass, and the mass-loss 
rate may vary, becoming increasingly higher as the 
MS phase comes to and end. Also, as noted above, in 
several cases the bubble will not reach this radius but 
may stagnate at a smaller radius. Due to all these 
reasons the ratio given in O may be reduced by up to 
a couple of orders of magnitude. However clearly the 
mass swept up by the shell significantly dominates 
over the mass expelled in the wind, and the wind 
material is not significant for the dynamics. Due to 
the large mass of the swept-up shell compared to that 
of the wind, the subsequent evolution is generally 
contained within the main-sequence bubble. Thus, 
the radius of the bubble does not change significantly 
after the main-sequence stage, although the internal 
structure may undergo substantial change. 

The swept-up mass lies in a thin, dense shell sur- 
rounding the bubble cavity, which consists of an in- 
ner freely expanding wind region followed by a region 
of almost constant density. The total wind mass 
ejected over time t is Mt. The mass in the freely 
expanding wind is MR t /v w (Rt = radius of wind 
termination shock). Since v w t >> R t by the end 
of the MS phase, a lower limit to the average cav- 
ity density is obtained by assuming that the wind 
material is uniformly distributed: 

SMt _3 / 2M 2 / 3 p a \ 3/5 4/5 

Pbub ~ 4mRl h ~ Att 0.76 3 [ ) 

(6) 

which for the main sequence stage can be written as 

PbubMS = 1.5xl0- 28 M% 5 p% 5 t lc T 4/5 g cm" 1 

(7) 
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Although this is a lower limit, especially if the 
bubble stalls early in the MS phase, it can be seen 
that the density in the interior of the bubble is on the 
order of 10~ 4 — 10~ 3 particles cm -3 , orders of mag- 
nitude lower than that of the surrounding medium. 

The position of the wind-termination shock is 
of interest. An accurate calculation is provided in 
Weaver et al. 1977 and Chevalier & Imamura 1983. 
To obtain a simple approximate expression we as- 
sume that the shock is strong, the shock jump is a 
factor of 4, and the bubble density is approximately 
constant at the post-shock value. In order to ac- 
commodate the results of W77, which indicate that 
the density in the shocked wind increases by about 
a factor of 2 close to the contact discontinuity, we 
include a factor a, with 1 < a < 2. The total wind 
mass emitted over time t is Mt, and the mass of the 
freely expanding wind is MR t /v w . The difference 
between the two gives the total mass of the shocked 
wind region between R c and Rt, where it is assumed 
that R c ~ R s h- 



Mt- 
which gives 

Rr 



MR* 



47T , , 



p 3n 



4a \R t /t 



4aM 



-1/3 



(8) 



(9) 



Note that this equation is similar to equation 13 
in Chevalier & Imamura 1983, if we take M<x >> 1, 
and assume S2 = Rt/t. For our purposes this approx- 
imation is sufficient. We can simplify this further by 
noting that in general v w t/Rt » 1. Then we get: 



Rt 



AaR c 
3v w t 



1/2 



(10) 



For the parameters of the MS star, this gives 

D o AO 1/2 1/10 , 2/5 -3/10 „V3/10 

Rt = 2.43 a 1 v 2 ' 50Q t w '° p w ' MJ 7 pc (11) 

The small value of the wind termination shock 
shows that the hot shocked wind occupies most of 
the volume of the main sequence bubble. 

2.2. Post Main Sequence Phases 

Stars with M < 50M Q generally become a red 
supcrgiant (RSG) after leaving the MS. These stars 
have large envelopes and slow winds (10-20 km s _1 ) 
with a high mass loss rate of 10~ 5 to 10 -4 Mq yr _1 . 
This creates a high density region around the star, 



confined by the pressure of the main-sequence bub- 
ble. For a RSG lifetime of about 200,000 years the 
total energy input in the RSG phase is 



E 



RSG 



x 10 46 M_ 4 U20 *o.2 ergs 



(12) 



where M_4 is the mass loss rate in terms of 1O~ 4 M 
yr , V20 is velocity in units of 20 km s _1 , and time 
is in units of million years, to. 2 = 0.2 x 10 6 years. 

Note that the RSG stage does not generally form 
a wind-blown bubble, because the velocity of the 
RSG wind is much lower than that of the medium 
(MS wind) into which it is blowing. However it does 
lead to a new pressure equilibrium. Since the total 
ram pressure of the RSG wind is different from that 
of the MS wind, the position of the wind termination 
shock will change. The new position can be found by 
equilibrating the ram pressure (pv%,) to the thermal 
pressure in the MS bubble. According to W77 the 
thermal pressure in the MS bubble is given by 



P bub = 0.163 L 2 / 5 Pa ^ t" 4 / 5 



(13) 



Equating to the ram pressure of the RSG wind 
gives the position of the wind termination shock: 



R 



tRSG 



1/2 



MrsqVRSG 
47T P bub 

8.85 Mr1\ v R20 - 5 A 

„ " 2 / 5 n -3/10 
y MS 2 500 PWa 



(14) 



-1/5 



tMS% 5 PC (15) 



where in the last expression the subscript R refers to 
the RSG wind and the subscript MS to the main se- 
quence wind. It is not clear however if the RSG wind 
will always be able to expand out to the distance re- 
quired to attain the new pressure equilibrium, and 
therefore in some cases the pressure equilibrium may 
never be attained. 

The RSG wind with its low velocity expands a 
distance Rrsg, with wind density prsg 



Rrsg = k 4.2 V20 h.2 pc; 
Prsg = 2.81 x 10~ 23 M 



-4 W20 1 



(16) 
(17) 



where we have added a factor k > 1 to account for 
the fact that neither the transition, nor the change 
in velocity, from the MS to the RSG phase is instan- 
taneous. Given the size of the MS bubble, it is clear 
that the RSG wind region will generally be confined 
to a small fraction of the main-sequence bubble. The 
total mass lost during the RSG phase is 



Mrsg = 20 M-4 t 



0.2 



Mr. 



(18) 
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Thus although total energy of the outflow in the 
RSG stage is small compared to the MS stage and 
the subsequent Wolf-Rayet stage, a large amount of 
stellar mass may be lost in the RSG stage. 

Some stars may go from the RSG to a blue super- 
giant (BSG) stage, as in the case of the progenitor 
of SN 1987A. The bipolar structure seen around the 
object is often interpreted as resulting from the in- 
teraction of a BSG and RSG wind. This is discussed 
further in SJ31 The wind parameters in this case are 
uncertain, but the mass-loss rate appears to be lower 
than even the MS stage, while the wind velocity is 
intermediate between the MS and RSG stage. 

Solar metallicity stars above 30-35 M Q end their 
lives as Wolf-Rayet (W-R) stars, although they may 
go through a luminous blue variable (LBV) stage. 
The mass-loss decreases somewhat to about 10 -6 — 
10 _5 Mq yr _1 , while the wind velocity increases to 
2000 km s _1 . For a lifetime of 100,000 years, the 
total energy input in the W-R phase is then 

E W r = 4 x 10 49 M- 5 v% 000 i .i ergs (19) 

Although the total mass is less, due to the high 
mass-loss rate and wind velocity W-R winds may 
posses enough momentum to push out, and possibly 
break up, any dense shell surrounding the star from 
the previous intermediate wind stage, distributing its 
contents throughout the surrounding medium. Gen- 
erally they will have enough momentum to collide 
with the MS shell, sending a reflected shock back. A 
W-R wind termination shock will be formed where 
the thermal pressure of the shocked wind bubble 
equals the ram pressure of the freely flowing wind. 

The post-MS stages may add considerable mass 
to the bubble without increasing the volume much. 
However equations ([6]) and ([7]) imply that, even in- 
creasing the mass (Mt) by a factor of 30-40 results in 
number densities of order 10 -3 — 10~ 2 cm -3 . There- 
fore the density over the bubble interior is in general 
low for W-R bubbles. 

Multidimensional calculation (Garcia-Segura et 
al. 1996; Freyer et al. 2006; Dwarkadas 2006b) reveal 
the presence of hydrodynamic instabilities in many 
stages. In one calculation of the medium around a 
35 Mq star, Dwarkadas (2006b) found that both the 
RSG and subsequent W-R wind (expanding into the 
RSG wind) were Rayleigh- Taylor unstable (Fig la). 
These instabilities tend to break-up the RSG shell 
and distribute its material over the entire bubble. 
They may lead to the formation of blobs, clumps and 
filaments in the WBB. Due to density fluctuations in 
various stages, the bubble interior becomes turbulent 
by the end of the simulation. 



Other factors will considerably modify this sim- 
ple picture. Rotation of the star can lead to an 
increase in the mass-loss rate (Maeder & Meynet 
2000). Dwarkadas & Owocki (2002) showed that a 
star rotating close to its break-up velocity may emit 
a wind preferentially in the polar direction. Mass- 
loss rates, which depend roughly as the square root 
of the metallicity, may be reduced at lower mctallic- 
ities, although rotation can still lead to high rates. 
And the presence of a binary companion can signifi- 
cantly alter the evolution, as has been hypothesized 
for SN 1987A (Morris & Podsiadlowski 2005). 

3. SUPERNOVA EVOLUTION 

The above description clearly illustrates a few 
salient points: 

• The medium around most massive stars con- 
sists of a low density cavity created in the main- 
sequence stage, surrounded by a dense shell. 

• If the progenitor is a RSG then the density near 
the star may be quite high, as expected for a 
RSG wind density. If the RSG stage is followed 
by a Wolf-Rayet stage then this can distribute 
the RSG material over the entire WBB, and the 
density will be much lower. 

• In cither case the region close in to the star is 
usually a freely flowing wind, with a wind den- 
sity that generally decreases as r -2 . However if 
the wind parameters change substantially with 
time then the density dependence may change. 
Several authors have suggested that the den- 
sity of the medium surrounding SN 1993 J de- 
creases only as r -15 (Suzuki & Nomoto 1995; 
Mioduszewski et al. 2001; Bartel et al. 2002). 

When the star explodes as a SN the resulting 
shock wave will first expand in the surrounding wind 
medium, and then in the low density cavity. If the 
surrounding wind medium is a RSG then the shock 
wave will expand in a higher density medium, and 
therefore its luminosity, due to circumstellar inter- 
action (Chevalier & Fransson 1994) will be high. If 
the shock wave results from the death of a W-R star 
then the surrounding medium density will be much 
lower, and the resulting luminosity will also be lower. 

The shock wave will expand in the wind until it 
reaches the wind termination shock, after which it 
will continue to expand in the low density shocked- 
wind medium. If it was earlier expanding in a RSG 
wind then its luminosity will be expected to drop. 
If it was expanding in a Wolf-Rayet wind then the 
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luminosity will not change much, or may in fact in- 
crease slightly as the shock wave crosses the wind 
termination shock. 

Since most of the mass is in the dense shell, 
the interaction of the shock wave with the dense 
shell controls the relevant dynamics, which depends 
mainly on one parameter A, the ratio of the mass 
of the shell to that of the ejected material. Sim- 
ulations for various values of A have been carried 
out by Tenorio-Tagle et al. (1990) and Dwarkadas 
(2005). They show that values of A > 1 are dy- 
namically important. As the value of A increases, 
the kinetic energy transmitted from the ejecta to 
the shell gradually increases. The collision with the 
shell results in a transmitted shock expanding into 
the shell, and a reflected shock moving back into the 
ejecta. It also results in an increase in the emission 
from the remnant, especially the optical, X-ray and 
radio emission. The reflected shock moves back to- 
wards the origin, thermalizing the ejecta on its way, 
faster than the original SN reverse shock would have. 
For smaller A the SN shock wave eventually 'forgets' 
about the existence of the shell, and the solution 
resembles what it would be in the absence of the 
shell. The density structure changes to reflect this. 
An increase in A results in higher velocity reflected 
shocks, while the transmitted shock slows down con- 
siderably. In extreme cases the transmitted shock 
may be trapped in the dense shell for a significant 
amount of time, and the shock wave may lose consid- 
erable energy to become a radiative shock wave. In 
such a case it can go from the free-expansion stage to 
the radiative stage without ever going through the 
Sedov or adiabatic phase. The remnant is also con- 
fined to the shell, whose size as we have seen was set 
in the MS stage. 

The evolution of SNe in wind blown bubbles thus 
differs considerably from the self-similar solutions so 
often used to describe SN evolution in general. The 
radius and velocity of the remnant do not evolve in a 
self-similar fashion but may vary considerably once 
the shock expands beyond the freely-flowing wind. 
Furthermore the expansion parameter of the rem- 
nant S (where R oc t s ) is continuously varying, as 
opposed to a self-similar case where the evolution is 
constant. This is illustrated in Dwarkadas (2005). 

In multi-dimensions (Dwarkadas 2006a, b) the in- 
terior of the nebula shows signs of turbulence, with 
significant density and pressure variations. The 
shock wave expanding in this medium no longer re- 
mains spherical, but becomes wrinkled due to the 
continuous interaction with the pressure and density 
inhomogeneities. It develops a corrugated structure 
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Fig. 1. (a) The evolution of the W-R wind into the RSG 
wind region for a 35 Mq star. Both winds are R-T unsta- 
ble, (b) The SN shock wave expanding into the density 
inhomogeneities in the interior can becomes wrinkled and 
lose its spherical shape (Dwarkadas 2006a, b). 

(see Fig lb) , and its interaction with the surrounding 
shell no longer occurs all at once but in a piecemeal 
fashion, with some parts of the shock colliding with 
the shell before others. As mentioned before, the 
collision results in an increase in the luminosity. In 
this case the luminosity of some parts of the shell 
will increase first, followed by those in other parts 
of the shell. A similar scenario appears to be taking 
place in the ring around SN 1987A. 

This brief description encapsulates the basic 
properties of the evolution of SNe in the medium 
around massive stars. Further details are given in 
Dwarkadas (2005, 2006a,b). In the rest of this paper 
we would like to concentrate on applying the results 
described above to SN 1987A. 

4. SN 1987A 

SN 1987A provides one of the best, and most 
spectacular, opportunities to witness the evolution 
of a SN shock wave within a wind-blown bubble. 
The three ring structure around SN 1987A has been 
interpreted as an hour-glass shape bubble formed 
by the interaction of a blue-supergiant (BSG) and 
a RSG wind (Blondin & Lundqvist 1993; Sugerman 
et al. 2005). The relatively small size of the equato- 
rial ring (possibly the waist of a wind-blown nebula) 
has ensured that we can see the shock-shell collision 
unfold in a short time span after the SN explosion. 
SN 1987A has allowed us to refine our theories of 
SNe evolving in the winds of massive stars, while 
providing confirmation of many existing ideas, and 
catalyzing many new ones. 

Chevalier & Dwarkadas (1995) showed that the 
slow expansion of the radio source, its large size, and 
the almost linear increase in radio and X-ray emis- 
sion could be understood if we assumed that the SN 
shock wave was interacting with an HII region inside 
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Fig. 2. The density profile just after the start of the simu- 
lation. The various regions (free expanding wind, almost 
constant density bubble, HII region and equatorial ring) 
are marked. 

of the inner ring. In the years since this model has 
been verified by Dick McCray and his collaborators 
(Michael et al. 1998). This model modifies the pic- 
ture outlined above by including a high-density HII 
region inside of the outer shell, in the low density 
wind-driven cavity. A depiction of the surrounding 
medium into which the SN explodes (in the equato- 
rial region) can be seen in Fig. 2. 

The morphology and dynamics of the ejecta in 
SN 1987A are inherently three-dimensional. The 
gradual appearance of bright spots around the equa- 
torial ring (Lawrence et al. 2000), and the changes 
in the X-ray and optical luminosity indicate an as- 
pherical shock interacting with several protrusions, 
formed perhaps by instabilities, emanating from the 
equatorial ring. All this presents a very complex 
morphology for the modeler. 

Yet there are several aspects of SN 1987A that 
can be understood by resorting to a simple, spher- 
ically symmetric model. This includes the radius 
and velocity of the expanding shock wave, the re- 
emergence of the X-ray and radio emission at around 
1100 days after explosion, and the gradual increase 
in the X-ray emission. We present here a toy model, 
a spherically symmetric calculation that captures the 
basic idea of Chevalier & Dwarkadas (1995), and 
compute from this the shock dynamics, the hard X- 
ray emission and radio emission, which we compare 
with observations. The hydrodynamic simulations 
have been carried out using the VH-1 code, a spher- 
ically symmetric multi-dimensional finite-difference 
hydrodynamic code. 

Figure 2 shows the initial conditions for the run. 
The SN ejecta, the wind termination shock of the 



BSG wind, the HII region and the equatorial ring 
are shown. The various parameters used for the run 
were: BSG wind mass-loss rate M = 5 x lO -9 M0/yr, 
BSG wind velocity = 550 km s , radius of wind ter- 
mination shock = 1.5 xlO 17 cm, inner radius of HII 
region = 4.5 xlO cm, inner radius of equatorial 
ring = 5.9 xlO 17 cm. The wind density is obtained 
from p w = M I ' (Aitr 2 v w ). The density is assumed to 
jump a factor of 4 at the wind termination shock, 
after which it remains constant up to the HII region, 
where the density jump is a factor of 100. The HII 
region density is assumed to increase steadily until 
the inner ring radius, whose density starts at 10~ 20 
g cm~ 3 and increases linearly. The SN ejecta den- 
sity varies as r -9 . These parameters were initially 
selected based on observational and theoretical work 
by Lundqvist (1998), Michael et al.(1998), Manch- 
ester et al. (2002), Park et al. (2006), and then mod- 
ified to obtain a better fit to the observations. 

In a future paper we will discuss the detailed hy- 
drodynamics of the evolution. Herein we concen- 
trate on the main results. The radius and velocity 
of the forward shock wave are shown in Fig 3. These 
quantities are in good agreement with 2 important 
observations: (1) Radio data from Manchester et 
al. (2002) which suggests that the shock velocity re- 
duces abruptly from a large value (> 35000 km s _1 ) 
to a value of around 3000 km s" 1 around day 1100; 
and (2) X-ray observations (Park et al. 2005, 2006) 
which indicate that around day 6000-6200 the shock 
velocity decreased to about 1500-1600 km s . 

Figure 4 shows an approximate computation of 
the radio emission from the remnant. We follow the 
prescription of Chevalier (1982), assuming that the 
radio emitting region lies between the forward and 
reverse shocks, with the optically thin radio lumi- 
nosity being given by 

L v oc 4nR 2 ARKB^/ 2 ^-^/ 2 (20) 

where L v is the radio luminosity at frequency u, 
AR is the thickness of the synchrotron emitting re- 
gion, B is the magnetic filed, and the distribution of 
accelerated particles is assumed to be a power-law 
N(E) = KE~~t. We further assume that both B 
and K scale simply as the hydrodynamic variables: 

U B = e B U th ; U re i = e r U th (21) 

where Ub is the magnetic energy density, U re i is the 
relativistic particle energy density, and Uth is the 
thermal energy density. The value of 7 is taken to 
be 2.7 (Manchester et al. 2002). 

In Fig 4 the solid line shows our prediction of 
the radio luminosity, while the symbols represent the 
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Fig. 3. The time evolution of the radius (top) and veloc- 
ity (bottom) of the SN forward shock wave. 



observed data at 1.4 GHz (Manchester et al. 2002). 
While the fit over the observed range is within a fac- 
tor of a few, we find it difficult to get the correct 
linear slope. In our calculations the emission is usu- 
ally increasing quadratically or faster, and by the 
time the shock hits the ring the radio emission be- 
gins to increase at an increasingly fast pace, which 
is not found in the observations. In our exploration 
of the parameter space we have not found a suitable 
set of parameters which generate a linear increase in 
radio luminosity. 

This model is unusual in that the best fit is ob- 
tained when either the magnetic field or the param- 
eter K representing relative particle energy density 
is assumed constant, with a slight preference for the 
latter. The reason why either quantity being con- 
stant works is due to the scaling in equation 1201 
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Fig. 4. The evolution of the radio luminosity with time, 
assuming that the parameter K is constant. 

The important parameter is the quantity KB^ 1+1 ^ 2 . 
With 7^3, this parameter becomes KB 2 . Coupled 
with the scaling in equation this implies that ei- 
ther K constant and B oc P 5 , or B constant and 
K oc P, both give KB 2 « Uth oc P, where P is 
the thermal pressure. Thus this behavior is simply 
a result of 7 being close to 3. It is not clear why a 
constant value for either quantity works best. 

Figure 5 shows a comparison of the hard X-ray 
luminosity computed using the CHIANTI code (Dere 
et al. 1997, Landi et al. 2006). According to Park et 
al. (2005,2006) the soft xray emission may be com- 
ing from the interaction of the shock wave with the 
protrusions in the ring. Furthermore line emission 
forms a strong component of the soft xrays. These 
computations are beyond the scope of this paper, 
and we have concentrated only on the hard xray 
emission. We also assume that, for the given shock 
speed, the electron temperature is much lower than 
the ion temperature by a factor of 50 (see Ghavamian 
et al. 2006). This provides a reasonable fit to the X- 
ray emission. The solid line in Fig. 5 is the emission 
calculated from our simulations, the data points are 
from Park et al. 2006. The comparison is reasonably 
good before the shock hits the ring, but there is a 
significant bump in the luminosity upon shock-ring 
interaction, something that is not seen in the data. 
Most of the hard X-ray emission appears to arise 
from the reverse-shocked ejecta. 

There are several shortcomings to this simple 
spherically symmetric approach, and the subsequent 
computation of the radio and X-ray emission. Rel- 
ativistic electrons may not be injected in a power- 
law fashion, the magnetic and/or relativistic par- 
ticle energy density may not vary as the hydrody- 
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Fig. 5. The variation of the hard X-rays with time. T e 
= 0.02 Tj, where 'e' = electrons and T = ions. The 
emission is mostly coming from the reverse-shock region. 

namic variables, or the variations may be different 
as the shock traverses the various density structures. 
The spectral index varies somewhat over the evo- 
lution (Manchester et al. 2002). The expansion is 
not spherically symmetric, so a spherically symmet- 
ric model may not necessarily work well. There is the 
possibility that with asphcrical expansion light travel 
time effects may come into play. The ratio of the 
electron to ion (post-shock) temperature, assumed 
constant here, may vary throughout the evolution. 
There will be contribution to the Xrays from the 
shock wave hitting the protrusions from the equa- 
torial ring, which will be discussed in a more com- 
prehensive model in future. We cannot hope to get 
much better agreement from such a simplified model, 
and it is satisfying that it does give reasonably good 
agreement with the observed data. 

In a future paper, we will discuss in more detail 
the hydrodynamics of the evolution and the reverse 
shock dynamics, which space constraints do not per- 
mit us to include herein. We will investigate more 
thoroughly the emission from different regions, as 
well as the soft xray emission. What is also needed 
is a direct computation of the medium around SN 
1987A, including the ionization from the star and 
the formation of the HII region, at least in the equa- 
torial plane. These will be addressed in future. 
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